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Abstract. A novel family of exactly solvable quantum systems on curved space is presented. 
The family is the quantum version of the classical Pcrlick family, which comprises all maxi- 
mally superintegrable 3-dimensional Hamiltonian systems with spherical symmetry. The 
high number of symmetries (both geometrical and dynamical) exhibited by the classical sys- 
tems has a counterpart in the accidental degeneracy in the spectrum of the quantum systems. 
This family of quantum problem is completely solved with the techniques of the SUSYQM 
(supcrsymmetric quantum mechanics). We also analyze in detail the ordering problem ari- 
sing in the quantization of the kinetic term of the classical Hamiltonian, stressing the link 
existing between two physically meaningful quantizations: the geometrical quantization and 
the position dependent mass quantization. 
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1 Introduction 

A Hamiltonian system of A degrees of freedom is said to be completely integrable, in the 
Arnold-Liouville sense, if it possesses A functionally independent globally defined and single 
valued integrals of motion in involution pp. In analogy to the classical mechanics, a quantum 
mechanical system with A degrees of freedom is called integrable if it is equipped with a set of 
A — 1 algebraically independent observables (self-adjoint operators densely defined on a suitable 
Hilbert space) all of them commuting with the Hamiltonian itself. The maximum number of 
independent observables commuting with the Hamiltonian is 2 A — 2; in that case out of this set 
one can extract at least two (complete) sets of A commuting observables, both containing the 
Hamiltonian. Such a system is called maximally superintegrable (M.S.). 

From a classical point of view these systems are characterized by the so called "closed or- 
bit property" for bounded motions, and from a quantum point of view these systems present 
a maximal degeneracy in the discrete spectrum. Over the years the M.S. systems have been ex- 
tremely useful models in the description of physical phenomena, as for instance the well known 
oscillator and Kepler radial problems. Recently, exactly solvable 1-dimensional systems have 
been proposed, above all in quantum mechanics, as models describing electronic properties in 

*This paper is a contribution to the Proceedings of the Workshop "Supersymmetric Quantum Me- 
chanics and Spectral Design" (July 18-30, 2010, Benasque, Spain). The full collection is available at 
|http://www.emis.de/journals/SIGMA/SUSYQM2010.html1 
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condensed matter physics, for instance in semiconductors and quantum dots [2j [U [5j O [7] ; the 
main feature of these systems is the presence of an effective mass in the Schrodinger equation, 
that can be considered as a position-dependent particle mass. Since, from a classical point of 
view, there exists a complete classification of all radial M.S. systems (flat case Bertrand theo- 
rem [8], curved space Perlick theorem [9]), the main goal of this paper is to "translate" these 
classical systems in the quantum context. Speaking about quantum M.S. systems many authors 
have recently stressed the existing relations of such quantum systems with systems studied in 
supersymmetric quantum mechanics (SUSYQM) [10] and in fact in this paper we will achieve 
our goal by using the powerful tools of the (SUSYQM) in particular the intertwining techniques 
that in the end allow us to obtain the spectrum and the eigenfunctions in a pure algebraic way. 



2 Perlick's family of classical M.S. systems 

Let us consider the standard radial Hamiltonian describing a particle subject to a central force: 



H(p r , p e , P+, r, e, <j>) = \p 2 t +^2+ v(r). 

A theorem proved by Bertrand in the nineteenth century [8] establishes that the only radial 
potentials such that all bounded orbits are closed turn out to be: 

V(r) = wV, V(r) = --. 

r 

We can obtain a wider class of such systems if we take in exam the more general Hamiltonian: 

H(P r , P e , P h r,6,<p) = T + V(r) = ^P 2 + ^ + V(r). (2.1) 

The kinetic term of this new class of Hamiltonian systems can be regarded as that of a particle 
moving according to a geodesic motion in a complete Riemannian manifold K 3 with the metric: 

fir) 

Considering this new scenario, Perlick managed to find a classification of the pair of functions 
(f(r);V(r)) admitting stable closed orbit property [9]. This classification is indeed the most 
general [TT] and it contains as particular cases other systems with the closed orbit property as 
for example the so called multifold Kepler systems and Darboux III system [12]. The Perlick's 
systems can be classified in two parametric families that can be viewed as a deformation of the 
two original Bertrand systems: 

Family I (Kepler type) 



iir 2 /I 
ds 2 = p2{l + k2r2) + r 2 (d9 2 + sin 2 (0)# 2 ) , V(r) = - + k 2 + G (2.2) 

Family II (oscillator type) 

9 2(1 - Dr 2 ± J(l - Dr 2 ) 2 - k 2 r 4 ) , 9/ , „ x 9N 
dS = mi-Dr 2 ) 2 -^) ^ + r 2 (d9 2 + sin 2 (9W 2 ), 



V(r) = GT 



1-Dr 2 ± V(l - Dr 2 ) 2 - k 2 r 4 ' 



In the above formulas, G is just an additive constant, and indeed it could be harmlessly removed. 
However we prefer to keep it in order to be consistent with the formalism used in our previous 
papers |13| 111]. These systems were proven to be M.S. in a rigorous way in 2009 |13j and define 
a complete classification of all radial problems in a Riemannian manifold with radial symmetry. 
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3 The generalized Kepler problem 

As said before, the Hamiltonians with radial symmetry turn out to be of physical interest above 
all in the quantum context. So, this section will be devoted to give a short introduction of the 
quantization problem of systems defined in a non flat space and to provide the tools needed to 
solve them exactly. The main problem is how to define the quantum analog of the kinetic term 
in (|2.1|) . since there is an order ambiguity in its quantization T(r,p) —> T(r,p). This problem 
always arises when in the kinetic energy term there is a dependence on the coordinates |144 [31 [6] , 
which comes from the nonzero curvature of the space. So, our Hamiltonian H = fP 2 + V(r) 
could be quantized for example as: 

-f a d r f l - a - b d r f b + fcil + V{T ), a, b e R. (3.1) 

Since our goal is to preserve the M.S. property of the classical Per lick's systems at the quantum 
level, we will set aside, for a moment, the choice of the parameters a, b and start directly from 
a M.S. quantum system. To this aim we choose the Perlick I space with the parameter j3 = 1; this 
system is universally known as the generalized Kepler problem, namely the Kepler problem on 
a space of constant curvature [T5J [161 H3 UHl [191 EQ] ■ As many other M.S. quantum systems [TO] , 
this system is indeed linked to a so-called shape-invariant system. 



3.1 Shape-invariant systems and intertwining technique 

Let us introduce the 1-dimensional hyperbolic Kepler Hamiltonian, namely the radial part of 
the 3-dimensional Schrodinger equation on a space of constant negative curvature ("hyperbolic 
space" [32]): 



2 & h(q ~ 1) o cosh(fc 



H q = -af + . / - 2nk 



r 



q " r ' sinh 2 (/cr) sinh(/cr) ' 

It should be noted that the parameter k is linked to the scalar curvature R through the relation 
R = —6k 2 and in fact in the limit k — > we recover the radial part of the usual flat Kepler 
problem. This Hamiltonian system is a member of a small, though extremely relevant, family 
of systems called shape-invariant systems: the reason for this name will be clear soon. 

These systems can be solved easily by the "intertwining technique" |21j . which can be seen 
as a generalization of the Dirac ladder-operator technique. To clarify this point, let us define 
the "prepotential" function: 

WJr) = -—r + qlnismhikr)), q 2 < ^, q > 
q k 

and the following two ladder operators: 

A\ = d r + W' q (r) = d r -^+ qkcoth(kr), 

-A q = d r - W' q {r) = d r + - - qk coth(fcr). 

The Hamiltonian turns out to be factorized by the ladder operators: 

6\ A -ir _ f e -_tL_h 2 n 2 

■"■q-n-q — 12 q fc <J) t q — ^2 " ' 

Since A q e Wq ( r > = we have: 

A\A q e w ^ = (H q - e q )e w ^) = Hq e w ^) = €q e w ^\ 

therefore we can read e^M = ^ ( r ) 

as an eigenfunctions of the Hamiltonian H q . 
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The crucial property turns up by intertwining the two ladder operators: 

A q A q = Hq+i — € q , 

which shows that we get the same Hamiltonian system, apart from the value of the parameter q 
that turns out to be changed to q + 1. This is indeed the exceptional property of the shape 
invariant systems. This property allows us to find in a purely algebraic way the spectrum and 
the eigenvectors. In fact, let us define the following relations: 

Al 1+1 A q+ iip q+1 (r) = (H q+l - e q+1 )ip q+ i(r) 
= {A q A\ + e q -e q+l )^ q+1 = 0. 

Multiplying by A J from the left: 

(A\A q + e q -e q+1 )(Ali> q+1 ) = 
and finally: 

H q {A% +1 ) = e q+1 {A\% +l ) 

therefore ipi, q (r) = A q ip q +i can be considered the first excited state of H q . 

By iterating this procedure it is possible to gain the whole set of bound states eigenfunctions 
as well as their eigenvalues. In particular the eigenvalues read: 

H q ip n ,q(r) = e q+n ip n)q {r) = - f ^q^p + fc2 (<? + n ) 2 ^j ^n,q{r), n £ N, 



while, up to a normalization factor, the finite set of the eigenfunctions will be: 

n-l 



(r), n < AW, (q + N max y 



i=0 



Ik 



AS, 



As stated at the beginning of this section this system can be seen as the k deformation of the 
radial part of the hydrogen atom Hamiltonian: in fact, in the limit k — > 0, its eigenvalues (and 
eigenfunctions) turn out to be exactly those of the hydrogen atom: 

e s = 7T, s = q + n. 

s z 

3.2 Hyperbolic Kepler and Perlick systems 

The radial part of the hyperbolic Kepler system is an exactly solvable 1-dimensional quantum 
system. In this subsection we are going to see how to link this 1-dimensional problem to one 
of the 3-dimensional problems associated to the metric (|2.2p . This link can be found easily by 
making the following point transformation: 

, sinh(/cr) 
r = 



yielding 



l 2 smli \kr) 



2 V ' r 2 ' 1r 

namely the quantization (|3.ip with parameters a = i, b = 0. 



r 
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The scalar product will change to the new one: 

dr' 



{lpn,q\lpn',q>) = / ( r ) ^n',q' (r)dr / £ gt^i',?'^' 



Vl + k 2 r' 2 



Finally, in order to get a geometrically coherent Hamiltonian system, we perform a similarity 
transformation that (of course) doesn't change the spectrum, but modifies the expressions of 
the eigenfunctions 

r' 2 dr 1>*(f) (I - A ^ n Jr') 



lo Vl + k 2 r' 2 r ' 
We define this new operator as the Laplace-Beltrami (LB) Hamiltonian: 

H LB = ±Hr> = -\ ((l + k 2 r> 2 )d 2 + 2 -d r , + 3k 2 r%^ + - ,^7k~ 2 - % 

with eigenfunctions: 

1 n— 1 

^( / )=^n^( / )^ + " (r,) 

orthogonal to each other with the new weight function w(r') in the semiline r' £ [0,oo): 

u;(r Jar = — =. 

VI + £;V 2 

The problem solved above is still a 1-dimensional problem on the semiline. Then in order 
to establish an exact correspondence between this 1-dimensional problem and the generalized 
Kepler problem we have to "upgrade" this 1-dimensional problem to a 3-dimensional one. We 
define the standard 3-dimensional angular momentum operator in angular variables: 

L 2 = -(d 2 + cot(8)d e + -±—d; 
V sin (p) 

and their spherical harmonic eigenfunctions: 

L 2 Y^ m (Q^) =l(l + l)Y l , m (9,ct>). 
It is possible now write down the 3-dimensional Hamiltonian operator: 

--2 n k 2 



#LB3d = ~\ ((1 + k 2 r' 2 )d 2 + ^d r , + SfcVfl^ + ^ - /^f + fc 2 - y , 

H^ nHl {r')Y l)m {e,<P) = (- f - k 2 (n + l + l) 2 \ ^ (/^ 

\ 2(n + l + iy 2 / 

We can recognize this Hamiltonian to be exactly the quantum version of the 3-dimensional gene- 
ralized Kepler system in which the kinetic part turns out to be the Laplace-Beltrami operator 
of the Perlick metric of type I with j3 = 1. This solves in a natural way the ordering problem of 
the quantization 



1 ■ • fl k 2 

H L B3d = -^^d i {g l ' J ^gd j )-fi^ : ^ + k 2 -—, g = detg itj . (3.2) 

We underline also the presence of a (maximal) degeneracy in the spectrum, a property which is 
characteristic of a M.S. quantum system 

H-LB2,d^n,l+l{r')Y^ m {9,(j)) = H LB 3di>n>,l>+l(r')Y l)m (6,(j)), 

n' / n, I' ^ I, n' + 1' = n + I. 
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4 The conformal Perlick problem 

Let us rewrite the general Perlick I metric in a conformal way performing a point transformation: 
dr 2 2 

ds2 = r 2 {r ^ A _ k 2 rP) 2 ^ dr2 + r2d ^)- 

The classical Hamiltonian related to this new system of coordinates turns out to be: 

r 2 (r-^-fcV)V „ 2 , L 2 
^ = § + 72" 

Seen from this perspective, the physical meaning of this classical system appears to be twofold: 
in fact we could regard it as describing either a particle with a constant mass in a curved space 
or a particle with a non-constant mass m = r2 ^ r -$ 4 _ k 2 r $^2 i n a flat space. 

In Section [3] we have investigated the hyperbolic Kepler system, namely the particular case 
of the Perlick I systems with (3=1. Therefore in this case the transformation (|4.1|) is: 

. 2f 

r - 



1 — k 2 f 2 

In this new coordinate system the Hamiltonian (|3.2p turns out to be: 

-2~ o „(„ 1\\ £.2 

T 



,2~2^2/^2 2fc 2 f _ 2. g(g-l)\ (\ k 2 r\ k 2 

H^ = --(i-k 2 r 2 ) [d 2 +T ^ df+ - df -m_iy^_ + _y_ 



The above Hamiltonian is of course the LB quantization of the generalized Kepler in the new 
coordinate system. Now the question is: it is possible to get a "physical", though "nongeometri- 
cal" quantization of this Perlick I system without loosing the M.S. property? The answer to this 
question is affirmative, and indeed in the following we will introduce another M.S. quantization 
of this problem. 

To this end, let us first rewrite the scalar product in this new variable: 

, f 8r 2 



(lpnMn',1') = / ; = K ^*n,lir')^n',l'{r')dr' / - T^^^n^f)^^ {f)df 

Jo \/l + k 2 r' Jd \ l - K r ) 



with two possible domains: 
D = r e 



i ] ; (Loo ) . 



kj 7 \k' 

Following the strategy of the previous section we get: 



we define now the new Hamiltonian: 



II 



vm 



l==H hB ^l-k 2 r 2 , (4.2) 



\/l — k 2 f 2 

Vi i2-2\lf&, 2 * q(q-l) \ f 1 k 2 r\ k 2 
H vm = --{l-kr) ( d f + -„dr - J - li f - f + — ) - 
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We recognize easily the differential part of this operator as the radial part of the Laplacian 
operator; this is the so called Schrodinger quantization |22j which amounts to solve the order 
ambiguity in the quantization of the classical system putting the conformal factor in front of the 
differential part, in other words the quantizations (|3.ip with parameters a = 1, b = 0. Therefore 
in this case the Laplace-Beltrami quantization and the Schrodinger quantization turn out to be 
linked through a similarity transformation, up to a constant shift ^- proportional to the scalar 
curvature of this system. The set of eigenfunctions and eigenvectors are linked to each other, 
through easy algebraic manipulations. In Section [3] we showed the way to get by an algebraic 
computation the whole set of eigenfunctions of this Hamiltonian operator. 

Before closing this section we wish to summarize all the algebraic manipulations we have 
carried out, and to display the general eigenfunctions written in a smarter way in term of Jacobi 
orthogonal polynomials P(n,a,p){ x ) 



2/jtanh 1 (fcr) 



fc („+i+i) r n+l / ^ _ k ( n + l + -Q2 ^ _ k ( n + I + ]_)2 \ / 1 + k 2„2 



n 



(1 _ k 2 r 2y+l+\ \ ' k{n + l + l) ' k(n + l + l) J\ 2kr 
Al , , / ^ 2 k 2 (n + l + l) 2 k 2 \ ( . . 

H ^ r) = ~ U + / + + 2 " ¥ ) ^ (r) ' 

Remark 1. Laplace-Beltrami vs Schrodinger quantization 

The above link between the Laplace-Beltrami and the Schrodinger quantization is far from 
being a lucky occurrence. Let us introduce a general conformal metric: 

ds 2 = f(xi,X2,xs) (dx\ + dx\ + dx 2 ^ . 

The scalar curvature of this conformal metric turns out to be: 

R(X 1 ,X2,X 3 ) =2^ ^3 • 

The classical Hamiltonian of a free particle in such space is: 

T = (P 2 + P 2 + P 2 ) 

Computing the Laplace-Beltrami quantization we get: 

while the Schrodinger quantization leads to: 

r vm = -^(v 2 ). 

Now by a direct computation it is possible to verify that the two quantum operators are related 
by a similarity transformation, apart for an additive function which bears a deep geometrical 
meaning, inasmuch as it represents the scalar curvature of the space: 

—f vm f* = Tlb + —R(xi,X2,X S ). 
fx 16 

Therefore the spectrum of the two different quantizations coincide up to a constant shift only in 
the case of constant curvature: otherwise the scalar curvature term becomes a radial potential 
term and cannot any more be considered as a shift in the eigenvalues. 
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5 Non-constant curvature quantum Perlick systems 

The machinery showed up to now plays a fundamental role in the quantization of the general 
Perlick problem I, namely the metric (|2.2p with /3 ^ 1. The general case /3 7^ 1 is, geometrically 
speaking, a highly nontrivial space because the scalar curvature is no more constant as it was 
in the generalized Kepler: 



R(r) = - \ ( (p 2 - 1) (k*r 2 ? + r- 2 ?) + 2k 2 (l + 5/3 2 ) ) . 



In the previous sections we performed the quantization following two different schemes, but 
since the scalar curvature turns out to be position-dependent the quantizations are no more 
similarity-equivalent and we need to make a choice. The most suitable choice, which keeps the 
solvability property, turns out to be the quantization (|3.1I) with parameters a = 1, b = 0: 

r 2^ + k 2 rP) 2( 2 a L 2 \ 11, _ ~ 

H ™ = '—j^ L U + ~dr ~ ^ j ~ § (r * + fcV). 

Performing the standard separation of variables: 

^nlm(r,9t,9 2 ) = <t>{r) n ,iY l>m {e u 6 2 ), L 2 Y l , m (9 1 ,9 2 ) = l(l + l)Y l , m (9 u 8 2 ), 
the Schrodinger problem is reduced to: 

- r2(r "y r ' )2 (d 2 + \ - ^11) - £(rH» + *V)) M r) = E(n, l^r). 

Now let us apply some algebraic manipulations already used in the past sections: 

1 

r = r , a=-, 



(l-fcV 2 ) 2 /«2 l + a fl a 2 l(l + l)\ 11 (\ 



.0—1 - 



;i - k 2 r' 2 ) 2 f a2 2„ a 2 /(i + l)-i^i\ f,/l 



# = ^A^- = - + " ,^ 4 j - f (p + * V J ' 

This turns out to be exactly the Hamiltonian operator (|4.2p if we set: q = al + . Performing 
this substitution we found the new spectrum: 

E(„,,) = 'J? _ ^ + al + ^ ? e 

K ' 2{n + al + ^) 2 2 T 8 

Let us consider a as a rational number, i.e. a = mi,m 2 € N. The two different wave 
functions Vv> Vv,Z' with n' = n — sm 2 , I' = 1 + sm 2 (s E N) are associated with the same energy 
eigenvalue, and therefore the degeneracy is still present in the general case /3 7^ 1 (but rational). 

Now it is crucial to underline that, since n, I are natural numbers, the degeneracy of the 
spectrum can be preserved only if a is a rational number too; but this is exactly the restriction 
requested also in the classical case for the closed orbit condition. 

Of course it is still possible to translate these results to the Laplace-Beltrami Hamiltonian by 
a suitable similarity transformation. Clearly, in the LB case we would find a slightly different 
Hamiltonian with a potential correction term proportional to the scalar curvature: 

jH vm f = H LB + cR(r). 
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6 Conclusions 

In the present paper, the first of a series on the radial quantum M.S., we have completely solved 
the first family of the quantum Per lick's systems with parameter k 2 > 0, giving contextually 
a few insights about the ordering problem. The case with k 2 < is also doable and pretty much 
similar to the one investigated in the present paper; the main difference consists just in the set 
of bounded states that turns out to be infinite. 

All the problems related to this family turn out to be solvable by the factorization technique 
(SUSY quantum mechanics). We have also provided the general solutions of the eigenfunctions 
in terms of the Jacobi orthogonal polynomials. 

The property of solvability and degeneracy in the spectrum holds indeed for the second 
family as well, providing therefore a complete classifications of all the radial quantum M.S. 
Hamiltonians. 

The results about the second family will be presented in a forthcoming paper. 
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